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(57) Abstract 

A method of NQR testing a sample having at 
least one NQR property which varies with a given 
environmental parameter comprises applying excit- 
ation to the sample to excite NQR resonance and 
detecting the NQR response signal. The excitation is 
such as to excite an NQR response signal which is 
detectable during the detection step over a predeter- 
mined range of the environmental parameter. The 
invention also relates to a method or, and apparatus 
for, NQR testing for the presence of selected nuclei 
(particularly nuclei of integral spin quantum num- 
ber, such as ,4 N) in a sample. It also relates to a 
method of determining the spatial distribution of 
temperature within a sample. 
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METHOD OF AND APPARATUS FOR NOR TESTING 

This application claims priority from United Kingdom Patent 
Application No. 9106789.2, filed on 2nd April 1991, whose 
5 disclosure 1s Incorporated herein by reference- 

This invention relates to a method of, and apparatus for, 
Nuclear Quadrupole Resonance (NQR) testing a sample having at 
least one NQR property (particularly NQR resonance frequency or 
an NQR relaxation time, T x . T 2 , T 2e and T 2 *> which varies with a 
10 given environmental parameter, such as pressure, magnetic field 

or- more— par-t-1cular-l-y~temperature.— It~aUo— reiates_to._a__me.thod 

of, and apparatus for, NQR testing for the presence of selected 
nuclei (particularly nuclei of Integral spin quantum number, such 
as 14 N) In. a sample. It also relates to a method of determining 
.15 the spatial distribution of temperature within a sample. 

As an example, the Invention has application to the detection 
In the field of 14 N quadrupole resonance signals from the 
explosive RDX concealed In parcels or luggage or on the person, 
or deployed 1n explosive devices. As another example, it has 
20 application to the detection of concealed drugs, for instance at 
airports. 

NQR measurements have the advantage that they do not require 
the sample to be placed In a strong magnetic field, and therefore 
do not require the large, expensive and sample-size limiting 

25 magnet structures which are needed for nuclear magnetic resonance 
(NMR) measurements. 

Quadrupolar nuclei have a nuclear spin quantum number I equal 
to or greater than unity (I£-l). If they are half Integral, the 
quadrupole Interaction (in zero magnetic field) produces two 

30 doubly degenerate levels for the case of 1-3/2 and one allowed 
transition ()_-m/2>, three doubly degenerate levels for I-5/2 and 
two strongly-allowed (^3/2, 3/2*5/2) and one weakly-allowed 
(J4*b/2) transition, and so on. For Integral spin nuclei, the 
roost Important of which 1s X4 N with 1-1 , there are usually three 

35 levels and three transition frequencies, dropping to one for 
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nuclei In axial ly-symmetric environments. All there transitions 
have a characteristic frequency and relaxation time (or times) 
which can be used to Identify the substance under Investigation. 
These frequencies and relaxation times do not depend on the 
5 presence of other substances, provided that they do not have 
transitions In the same frequency range. 

Nuclear quadrupole resonance response signals are 
conventionally detected by means of pulsed radlofrequency (rf) 
radiation of the correct excitation frequency (v 0 > to excite the 
40— selected-tranmWn^tHrTesonan^ — a pu l se of ~ 
preset width t, rf field amplitude B lt and flip angle a generates 
a decaying signal Immediately following the pulse known as a free 
Induction decay (f.l.d.), and two or more pulses of pre-set 
widths and spadngs can generate echoes. The pulse width which 
15 produces the maximum f.l.d. for a given quadrupolar nucleus In a 
solid powder Is given the symbol t m and the corresponding flip 
angle the symbol (equivalent, for example, to a "90*" pulse in 
NMR). 

Usually in NQR testing a repetitive pulse train is used with 
20 spadngs between pulses, x, which depend on the relaxation times 
T lt T 2 and T 2e , and a number of f.l.d. 's and/or echoes Is 
accumulated to provide the required sensitivity of detection. By 
"sensitivity of detection" is meant the number of selected nuclei 
(or amount of material) that can be detected at a given level of 
25 probability in a fixed volume of sample over a given test time. 
This sensitivity may be further improved by Fourier 
transformation of the accumulated signals to give an absorption 
spectrum, the area under which is then measured by integration 
between appropriate limits. A paper by N.E. Ainbinder si &L 
30 (Advances in Nuclear Quadrupole Resonance, 2, 1978, 67-130) 
provides information on this background state of the art. 

The present invention seeks to provide an improved method of 
and apparatus for NQR testing. 

According to the present invention, there Is provided a 
35 method of NQR testing a sample having at least one NQR property 
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which varies with a given environmental parameter, comprising 
applying excitation to the sample to excite NQR resonance and 
detecting the NQR response signal, the excitation being such as 
to excite an NQR response signal which is detectable during the 
5 detection step over a predetermined range of the environmental 
parameter. 

The Invention Is based on our discovery that the NQR 
properties of a sample, such as resonance frequency (or 
frequencies), spin-lattice relaxation time (T x ), spin-spin 

10 relaxation time (T 2 >, effective spin-spin relaxation time (T 2e ) 
or free — Induction— decay — relaxation — time— (T 2 ^)< — may — vary- 
considerably with a given environmental parameter, such as 
temperature, pressure or magnetic field, and that 1t Is therefore 
necessary to take account of or compensate for these variations 

15 when carrying out NQR tests to cover a predetermined range of the 
environmental parameter. Such compensation can have the 
unexpected advantage of considerably more sensitive tests than 
have hitherto been achieved. 

Since the NQR testing would usually be carried out In the 

20 field, (such as at an airport), the given predetermined range of 
the environmental parameter may typically be the range of that 
parameter likely to be encountered in the field (possibly the 
national or global average range or maximum range of that 
parameter). If the parameter is temperature, the predetermined 

25 range may be ± 10 # C (possibly between 5'C and 25 # C) or ± 20°C 
(possibly between -10°C and +30°C). The range may be as much as 
from -30°C (corresponding to arctic conditions) to +40°C 
(corresponding to desert conditions). In some Industrial 
applications, considerably wider temperature ranges may be 

30 encountered. If the environmental parameter Is pressure, the 
predetermined range may be, for example, ± M (corresponding to a 
typical average daily pressure range) or ± 5X (corresponding to a 
maximum range). 

As used herein, the term "detectable" preferably connotes a 
35 signal/noise ratio which, within the . measurement time of the 
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detection step, Is significantly above the actual or expected 
■ n noise levels, the degree of significance being determined, for 
example, by standard statistical methods, such as the upper 
Student's t-distrl button. . . 

5 Take, as an example, the detection of NQR signals from the 
explosive RDX concealed In parcels or luggage or on the person, 
with the environmental parameter being temperature. The NQR 
- property resonance frequency 1s considered first. The three NQR 
resonance frequencies, v x ,v y ,v 2 , are all triplets due to crystal 
—-10 effects;— near^98-K—the-h-1ghTsr^ef-(^ ) occurs at 5047 (N A ), 
5192 <N 2 >, and 5240 (N,) kHz. Between 240 and 340 K, their 
frequencies in kHz follow closely the equations 
N i v Q (i) - 5148 -0.223T -0.000395T 2 

N 2 v Q (2) = 5277 -0.1 39T -0.000506T 2 (i) 
15 N 3 v Q (s) « 5332 -0.108T -0.000670T 2 

where T Is temperature. 

If the RDX were In airport luggage. It might be subjected to 
a temperature variation of -30 to +20°C about a mean of 20 C C 
(i.e. a variation from -10 to +40*C). Different samples might be 
20 at different temperatures, and/or the same sample might have a 
non-uniform temperature. Equation (1) predicts that vq (i) will 
vary significantly (between 5041 and 5062 kHz) for a temperature 
range of -10 to +40'C. The present invention can, for example, 
provide that the excitation Is such as to excite detectable NQR 
25 resonance at resonance frequencies between 5041 and 5062 kHz, 
that is, over a temperature range likely to be experienced with 
airport luggage. 

Considering now the NQR property spin-lattice relaxation time 
<T a ), we have discovered that the effects on the spin-lattice 
30 relaxation time (T a ) of temperature (and also pressure) 
variations can be very significant. For example, with the RDX 
sample described above, over the temperature range -10 to +40°C 
it Is believed that T A follows the expression 

T a - aT 2 + be c/kT (3) 
35 In which a and b vary according to the transition and c Is about 
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70 kJ rool" 1 . For H x 1n RDX T a varies, surprisingly, by a factor 
of approximately 8 between 5 and 25*C. 

He have further discovered that one particularly Important 
effect of the variation of T A Is on the ratio x/T lt where x is 
5 the pulse spacing between consecutive excitation pulse 
repetitions (sometimes known as the pulse repetition time). As 
an example, reference 1s made to Figure 1, In which the variation 
of signal, strength with flip angle (a) 1s shown for different 
values of x/T^ Figure 1 1s derived from equations derived by 

10 Vega (0. Chem. ' Phys., 1974, £1, 1099, Eq. IV-29) for sp1n-l 
— r- — nucle1-wh1ch-predict-th^ 

In a powder as a function of x for a given T 2 ; at a given ratio 
of x/Tj <<5) the signal Is predicted to pass through a maximum as 
the flip angle 1s varied, this maximum moving to flip angles less 

15 than and to lower signal strengths as the ratio of x/T 1 drops 
below 5, as shown In Figure 1. Suppose, at x/Tj-5, the maximum 
f.l.d., corresponding to a flip angle o^, has a relative strength 
of unity; at x/T^O.l, it has shifted to 0.3^ and has a 
relative signal strength of 0.25. For constant x and at a 

20 constant flip angle c^, a variation In x/J 1 between 1 and 0.1 (an 
Increase In l x by a factor of ten, corresponding, for N x in RDX, 
to a temperature Increase near ambient of roughly 20*C) gives 
rise to an approximately 70% loss in response signal strength. 
Such a loss may render the response signal undetectable if it is 

25 overwhelmed by noise. The present invention can accommodate the 
variation of T A with an environmental parameter such as 
temperature. 

Further, we have also discovered that the relaxation times 
T 2 , T 2e and T 2 * may vary with temperature or pressure, although 

30 in general It Is believed that this variation Is less marked than 
with T 1 . As 1s well-known, T 2 and T 2e are spin-spin relaxation 
times which are relevant to the generation of echoes. T 2 * 1s 
relevant to the generation of f.l.d. 1 s. One important effect of 
the variation of the spin-spin relaxation times 1s on the ratios 

35 x e /T 2 , x e /T 2e and x f /T 2 *, where x e is the pulse spacing between 
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the pulses used to generate one or more echoes and x f 1s the 
pulse spacing between the pulses (possibly at different 
frequencies) used to generate f.l.d.'s. . Variation of these 
ratios 1s Important because response signal strength 1s. a 
5 function of x e /T 2 , x e /T 2e and x f /T 2 * t an exponential being one 
form of function which is often encountered. Again., the present 
Invention can accommodate the variation of such relaxation times 
with an environmental parameter such as temperature. 

In putting the Invention Into effect, if the NQR property 
-^lO^under^consTderatlo^^ 

at least one excitation pulse (e.g. a simple or composite puise) 
at a single excitation frequency, the or each pulse having a 
power spectrum sufficient to excite a detectable NQR resonance 
response signal over the range of resonance frequency 
15 corresponding to the range of the environmental parameter. Thus 
all frequencies to be encountered lie within (preferably well 
within) the power spectrum of the excitation pulse or pulses. A 
disadvantage of this approach is the high rf powers required. 

Therefore, 1t is preferred that the excitation is applied at 
20 a plurality of different excitation frequencies so that a 
detectable NQR resonance response signal can be excited over the 
range of resonance frequency corresponding to the range of the 
environmental parameter. This can be achieved, for instance, by 
applying separate excitation pulses at different frequencies, or 
25 by frequency modulating the fundamental frequency of one or more 
excitation pulses, or by varying the frequency of each individual 
pulse. 

Now the signal strength for any transition varies according 
to the relative values of the off-set frequency Av 0 and the 
30 precessional frequency v 2 in the rf B x field:- 

35 2ir 

where y n is the gyromagnetic ratio (1 .932 x 10 7 rad s" 1 T 1 for 

14 N nuclei). The strength falls off significantly as soon as Av 0 

becomes much larger than v x . For the N A resonance at 5047 kHz in 
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RDX, the variation 1s shown In Figure 2 for a 60ps rf pulse of B x 
amplitude corresponding to a value of v x of about 9 kHz; the 
diminution 1n signal strength with increasing Av 0 becomes 
Important when Av 0 Is greater than roughly 9 kHz. 
5 In the present Invention, by applying the excitation to the 

sample at a plurality of different frequencies so that a 
detectable NQR resonance response signal can be excited over the 
range of resonance frequency corresponding to the range of the 
environmental parameter, the maximum frequency off-set can be 

10 reduced so- as to be kept within acceptable limits. From the 

— — point-of— view-of— keeping~the-ma^(-1mum~frfequency— off-set— wUhin 
acceptable limits, the number of different frequencies at which 
the excitation 1s applied should be as large as possible. On the 
other hand, a conflicting requirement 1s that the greater the 

15 number of frequencies, the longer may be the NQR test. 

In putting the Invention' Into effect, 1f the NQR property 
under consideration Is a relaxation time, then the excitation 1s 
preferably applied as a series of excitation pulses, at- a 
plurality of different pulse spadngs and/or so as to produce a 

20 plurality of different flip angles, so that a detectable NQR 
resonance response signal can be excited over the range of the 
relaxation time corresponding to the range of the environmental 
parameter. 

This Is based on our .discoveries, as described above, of the 
25 effects on the spin-lattice, sp1n-sp1n and f.1.d. relaxation 
times of, for example, temperature and pressure variations and of 
the effect that these have on signal strength. 

If the relaxation time under consideration Is T 1 , the pulse 
spadngs (In this case the pulse repetition times, t) and/or flip 
30 angles can be chosen according to the principles explained 1n 
relation to Figure 1 to excite a detectable NQR response signal 
over the entire T 2 range of Interest. If the relaxation time 1s 
T T 2e or T 2 *, since signal strength 1s a function of the ratio 
of pulse spacing to the relevant relaxation time as well as a 
35 function of flip angle, pulse spacing (In this case the spacing 
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referred to as t e or Tf ) may be varied In direct proportion to 
the relevant relaxation time, or else flip angle may be 
appropriately adjusted, to excite a detectable NQR response 
signal over the entire range of interest of the relevant 
relaxation time. 

It 1s advantageous that the NQR testing takes Into account 
variations both In resonance frequency and J, caused by 
variations in the environmental parameter or parameters. . To 
achieve this, preferably the excitation Is applied as a 
repetitWe-sertes-of-e*^^ a plurality of 

different excitation frequencies, the pulse repetition time being 
the same for each series, the flip angle being different for each 
'.. series. One particular advantage of this arrangement Is that the 
pulses associated with the different frequencies can be 
15 interleaved without a pulse at one frequency Interfering with the 
response signals generated by pulses at other frequencies. . 

If 1t is desired to take Into account variations in T 2 or 
T 2C , and the excitation Is applied as a series of excitation 
pulses at each of a plurality of different excitation 
frequencies, each series Including a plurality of pulses for 
generating echoes, then the spacing between each of the plurality - 
of the pulses (x e ) is preferably different for each series (and 
chosen according to the variation of T 2 or T 2e with the 
environmental parameter). The variation in Tj can thus still be 
25 compensated for by varying flip angle. 

If It Is desired to take Into account variations in T 2 *. and 
the excitation is applied as a series of excitation pulses at 
each of a plurality of different excitation frequencies, for 
generating detectable .f.i.d.'s. pulses at the different 
frequencies are preferably applied at staggered spaclngs (t f ) 
chosen according to the variation of T 2 * with the environmental 
parameter. x f may be some multiple of (say two or three times) 
T 2 *. such as allows, a reasonable proportion of the f.l.d. 
produced by the preceding pulse to be detected before the 
35 succeeding pulse Is applied, without being so long as to 
unnecessarily slow down the test. 



20 



30 
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These techniques for taking Into account variations In T 2 , 
T 2e or T 2 * are particularly advantageous when combined with the 
aforementioned technique for taking Into account variations In 
resonance frequency and T a . In this event, flip angle adjustment 
5 can be reserved for compensating for 1 1 variations whilst pulse 
spacing adjustment can be reserved for compensating for T 2 , T 2e 
or T 2 * variations. This leads to experimental simplicity. 

Preferably, the excitation is such that the NQR response 
signal /noise ratio would be substantially equal at at least two 

10 differing values of the environmental parameter. It • Is 

-preferable— that— the— slgna-1 /no1se-~ratlo— would—be-^substantlaLl-y^ 

equal at as many differing values of the environmental parameter 
as" possible, and that It would be substantially constant between 
these values (that Is, so that It would be substantially constant 

15 over the range of the environmental parameter). "Substantially" 
is to be Interpreted In the context of the large variations In 
signal /noise ratio which may be occasioned by changes 1n an 
environmental parameter. Thus by "substantially equal" and 
"substantially constant" may be meant a minimum signal/noise 

20 level within 50*, 601 or preferably 751 of a maximum signal/noise 
level. 

A particular advantage of the NQR response signal/noise ratio 
being substantially equal at at least two differing values of the 
environmental parameter and preferably substantially, constant 

25 over the range of the environmental parameter Is that it 
facilitates accurate NQR testing for the presence of selected 
nuclei in a sample, such as would be required in the example of 
RDX detection above. The invention envisages that an alarm 
signal would be produced If a predetermined threshold of 

30 detection were exceeded. This threshold might typically be 
derived from a number of tests carried out on a variety of 
substances other than the substance of Interest 1n order to 
determine expected levels of background noise. The threshold 
would typically be set above these expected levels (preferably 

35 several standard deviations above these levels). By providing 
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that the signal/noise ratio remains substantially constant over 
the range of the environmental parameter, the effect of the 
threshold will also be substantially constant- In other words, 
the sensitivity of detection will be substantially constant. 
. 5 Substantial equality in the signal/noise ratio could be 

obtained by "post-processing" of the detected NQR response signal 
, : by some form of noise filter. However, It Is more preferably 
obtained by "pre-processing", by appropriate adjustment of the 
response signal strengths. It will be understood that the . 
^10 "noise" component of "the sTgnaT/nolse ratio need not Be 
specifically measured. It can usually be assumed that noise will 
be virtually invariant across the range of the. environmental 
parameter.. If this Is so, then maintaining the signal/noise 
ratio substantially equal reduces to a requirement to maintain 
15 the NQR response . signal strength substantially equal. If noise 
actually varies across the range, then a suitably .varying 
threshold might be appropriate. 

The invention envisages that the NQR response signal /noise 
ratio may be rendered substantially equal at the at least two 
20 differing values of the environmental parameter by compensation 
for resonance frequency-, T 1( T 2 , T 2e or T 2 * variations. Indeed, 
compensation may be for a combination of or even all of such 
variations. 

If the NQR resonance frequency of the sample varies with the 
25 given environmental parameter, and the excitation is applied at a 
plurality of different excitation frequencies, it is preferred 
that the number and spacing of these frequencies 1s such that the 
NQR response signal/noise ratio is substantial 1y constant over 
the range of the environmental parameter. In the example given 
30 above, if vg (l) lies between 5041 and 5062 kHz over a temperature 
range of -10 to +40T, this range of vg (l) can be covered by the 
use of two excitation pulse sequences, one at a radio frequency 
of 5048 kHz and the other at 5055 kHz; the maximum off-set at 
any temperature within the specified range is. Av Q = 7 . kHz, and, 
35 with this value, either one pulse train or the other will give 
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optimum or near-optimum sensitivity of detection for a pulse 
width of 60 ps and v a £ 9 kHz, with the NQR response signal 
strength remaining substantially constant. On the assumption 
that noise 1s constant across the frequency range of interest, 
5 the signal /noise ratio and hence the sensitivity of* detection 
also remains substantially constant. 

The spectrometer In this example 1s tuned to a mean frequency 
of 5051.5 kHz, and the probe 0 factor and receiver band width are 
selected so that signal responses at 5041 and 5062 kHz are 
10 amplified to the same extent. After separate Fourier 
"tni^fofm^ 

pulse trains, the absorption mode signals from each are then 
separately Integrated between frequency limits pre-set to ensure 
that the majority of the spectral line lies within the 

15 integration limits anywhere In the frequency range from 5041 
to 5062 kHz; 1n this example, Integration limits of (say) 2 to 
10 kHz will accept a signal from one or other of the excitation 
frequencies of 5048, 5055 kHz. The lower limit of 2 kHz 1s 
selected to reduce the effects of changes In the base-line 

20 response signal. The Integrals from both rf sequences are then 
separately monitored, or added to give the summed output. It 
will be understood that, 1n the example described, excitation 
pulses at different excitation frequencies may In some 
circumstances excite a single resonance frequency to a 

25 significant extent. The summed output of such pulses 1s not 1n 
this experiment nearly as great as 1f pulses were applied 
separately In isolation from each other, and their outputs were 
then summed, because the first pulse saturates subsequent 
response signals (although clearly this may not always be the 

30 case). 

Preferably, for simplicity, the excitation Is applied at n-1 
equally spaced frequencies, where n Is the nearest Integer 
satisfying the equation 

n5 Av 0 /Av off , 

35 where n£2, Av 0 Is half of the value of the resonance frequency 
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range corresponding to the range of the environmental parameter 
and Av off is approximately equal to the precesslonal frequency 
(v,) (and is therefore a measure of the maximum frequency off-set 
above which significant reductions In response signal/noise ratio 
would be observed for a given resonance frequency). 

The lowest frequency may therefore be v 0 -Av 0 /2 + Av off , the 
next v 0 -Av 0 /2 + 2Av off , and so on. in the example given above, the 
integration limits for the Fourier transformed spectrum can then 
be set between (say) 2 kHz and *v off . The Integrals from all 
pu se-sequences-may-theh-be monitored consecutively or summed to 
obtain the final output. The receiver band width and Q-factor of 
the rf coil would be selected to provide constant sensitivity 
over the frequency range from v 0 -Av 0 to v 0+ *v 0 . 

If a relaxation time of the sample varies with the given 
15. environmental parameter, and the excitation 1s applied as a 
series of excitation pulses at a plurality of different pulse 
spaclngs and/or so as to produce a plurality of different flip 
angles, the pulse spaclngs and/or flip angles are preferably such 
that the NOR response signal/noise ratio Is substantially 
constant over the range of the environmental parameter. This can 
facilitate accurate NQR testing f or the presence of. selected 
nuclei in a sample, as discussed previously. 

If the NOR resonance frequency of the sample also varies with 
the environmental parameter and the excitation pulses are applied 
at a plurality of different excitation frequencies, the pulse 
spaclngs and/or flip angles associated with the excitation 
frequencies preferably are such as to produce an NQR response 
signal of substantially equal signal/noise ratio at the 
respective values of the environmental parameter at which the NQR 
resonance frequency equals the respective excitation frequency 
This affords a particularly simple way of compensating for the 
variation of resonance frequency and relaxation time. . 

As an example, the case where only two. excitation frequencies 
are utilised 1s considered. To ensure that the sensitivity of 
detection at two different frequencies with different spin- 
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lattice relaxation times remains the same, 1t 1s necessary to 
change the flip angle a according to Vega's equations, either by 
using different pulse shapes (especially widths) or different rf 
powers at the two frequencies. The method of selection Is 
5 Illustrated diagrammatlcally in Figure 1 for the steady-state 
signal generated by pulse trains of two different ratios of 
t/T 2 - If t/T x at one frequency and temperature 1s equal to 1, 
setting the flip angle to 0.17 a,,, will give the same response 
signal strength as at another frequency and temperature at which 

io -c/T x « 0,1 and the flip angle 1s set to 0.3 o^. ^Thus the same 

— -number-of-accumulati^ 

same Integrated signal strength and therefore generally the same 
signal /noise ratio. 

The flip angles selected will only be correct at the spot 

15 frequencies of the excitation pulses, but deviations at 
intermediate frequencies will be reduced In their effect If 
Integrated signal strengths from consecutive Irradiation 
frequencies are summed. In practice, Inhomogeneities 1n the. B 1 
field and varying flip angles caused by temperature shifts 

20 between samples and/or temperature variations within a sample 
will also ensure that the idealised predictions are not followed 
exactly. It may therefore be necessary In a given case to apply 
small adjustments to the predicted flip angles to allow for such 
effects with a given coil .geometry and NQR frequency. 

25 Whilst the excitation might suitably take the form of simple 

rectangular pulses, it may be preferable that the excitation 
includes at least one excitation pulse which is so shaped as to 
produce a response signal whose strength is substantially 
constant over a limited resonance frequency range. Simple 

30 rectangular pulses typically have a curved response signal 
characteristic like that shown in Figure 2; a relatively flat 
characteristic afforded by a suitably shaped pulse would have the 
advantage that the sensitivity of detection would be more uniform 
over the limited resonance frequency range. Suitable pulse 

35 shapes will be known to the skilled person. One particularly 
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suitable shape 1s an Hermlte pulse (see the paper by M. McCoy and 
W.S. Warren, J. Mag. Res., 1985, £5, 178). which has the 
additional advantage of requiring less power than a rectangular 
pulse to generate a given B x field. 

The invention extends to apparatus for NQR testing a sample 
having at least one NQR property which varies with . a given 
environmental parameter, comprising:- 

storage means 1n which 1s stored Information oh how the or 
each NQR property varies over a predetermined range of the 

TO env1ronmemil^aTSSet¥n ~ 

means for applying excitation to the sample to excite NQR 
resonance; 

means for detecting the NQR response signal; and 
means, responsive to the information stored in the storage 
15 means, for controlling the applying means to excite an NQR 
response signal which is detectable by the detection means over a 
predetermined range of the environmental parameter. 

The Invention further extends to a method of NQR testing for 
the presence of selected nuclei 1n a sample, comprising: 
20 deriving an NQR response signal from the sample in such a way 

that the signal/noise ratio would be substantially equal at at 
least two differing values of an environmental parameter; and 

producing an alarm signal from the NQR response signal in 
dependence on whether a predetermined threshold of detection has 
25 been exceeded. 

As alluded to earlier, by this arrangement, a particularly 
accurate and sensitive NQR test for the presence of selected 
nuclei can be carried out, since the signal/noise ratio will have 
the same significance In relation to the threshold value at the 
30 two or more differing values of the environmental parameter. 

This aspect of the Invention extends to apparatus for NQR 
testing for the presence of selected nuclei In the sample 
analogous to the method described Immediately above. All the 
features of the Invention as a whole described above apply to 
35 this aspect of the Invention. 
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The Invention extends to a method as aforesaid of NQR testing 
a sample having at least one NQR property which varies with a 
given environmental parameter, the method being for determining 
the spatial distribution of temperature within a sample, wherein 
5 a relaxation time <T lt T 2 , T 2C , T 2 *) varies with temperature, and 
such spatial distribution 1s determined according to the 
variation of relaxation time within the sample. 

This method 1s based on our discovery that relaxation time 
varies substantially with temperature. It can provide, a 

10 sensitive technique for temperature distribution Imaging. 

If - is preferFetT" tfiaf~ t"K¥~s"pa"t1'a"1 — var1at1on-d1str1but1on-of 

temperature 1s determined according to the variation of the 
spin-lattice relaxation time, T x , within the sample, since we 
have found that this relaxation time may be most sensitive to 

15 temperature. One preferred way of putting this method Into 
practice 1s to apply a series of excitation pulses to the sample 
In which the pulse repetition time Is varied whilst the flip 
angle 1s kept constant. The pulses would be applied, for 
Instance, 1n the presence of a weak magnetic field to provide 

20 spatial encoding of the data. Decoding of the data can be 
carried out by well-known techniques such as Fourier analysis. 
Tj data could be converted to temperature data via a look-up 
table. 

The spatial distribution of temperature may also be 
25 determined according to the variation of a spin-spin relaxation 
time within the sample. This method would be put Into practice 
by varying the spacing <t e > between pulses whilst keeping the 
flip angle constant, 1n order to encode spatially data relating 
to the relevant relaxation times. It is thought that this 
30 technique may apply to the times T 2 and T 2e but not to the time 
T 2 \ since In the case of T 2 * there is no relevant pulse spacing 
which can be varied. 

The invention extends to a method for detecting the presence 
within a larger article of a specific substance containing 
35 quadrupolar atomic nuclei, wherein the irradiation can allow for 
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resonance frequency shifts caused by *20«C temperature 
variations. It also extends to a method for detecting the 
presence within a larger article of a specific substance 
containing quadrupolar atomic nuclei, wherein the power spectrum 
5 of the Irradiation pulses provides substantial power within about 
0.1% of any frequency to which NQR resonance may be shifted by 
any environmental conditions likely to apply to the article. The 
Irradiation pulses may be at one or more frequencies equal or 
close to (maybe within 0.11 of) a resonant frequency of 

_ ^10_quadr-upolar-nuc-l*4-rw4th1n-the-siJbstaflce-to-be-detected- — r 

The substance may be drugs, for Instance heroin or cocaine, 
or explosives, for instance HMX, RDX, PETN, or TNT. The 
quadrupolar nuclei 1n these substances may be 14 N. The 
environmental conditions are such conditions as temperature, 
15 pressure or magnetic fields. One use for these methods Is In the 
examination of airline baggage or airfreight, where In a typical 
case ±10 kHz resonance frequency shifts may be caused by ±20 4 C 
temperature variations. For the detection of the explosive RDX a 
suitable resonance frequency is 5.191 MHz; rf drive signals. 
20 having a power spectrum of width 18 kHz at half height may be 
used to allow for the effects of the ±20°C temperature variation. 

For any aspect of the invention, apparatus features analogous 
to the described method features may be provided. 

Preferred features of the Invention will now be described, by 
25 way of example, with reference to the accompanying drawings, in 
which:- 

: Figure 1 1s a plot of the steady-state NQR response signal 
strength against flip angle for differing values of the ratio 
t/T x ; 

30 Figure 2 is a plot of NQR response signal strength against 

frequency off-set for a 60 ps excitation pulse (v x * 9 kHz); 

Figures 3 are timing diagrams for two embodiments of the 
invention; and 

Figure 4 is a block diagram of NQR apparatus used in the 
35 present invention. 
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Referring to the timing diagram of the first embodiment of 
the Invention shown 1n Figure 3(a). two radio frequency pulses 
(either simple or composite) of width t w <x) and t w <2> and of 
differing frequencies f x and f 2 are repeated at pulse repetition 
5 times x (x»t w ) which are much longer than the longer J x - say 
t-5T a - to ensure full signal recovery between pulses; phase 
shifts of alternate pulses or suitable combinations of pulses of 
width t w and 2t w may be used to eliminate probe ringing. After 
phase-sensitive detection and manipulation of the appropriate 
10 signals, the residual oscillations can be made to cancel and only 

the-true-NQR~response-s1gnal~1s-observed^ — Slnce-T-ls-muc-h-longer 

than either J x value, no l x compensation Is necessary (see Figure 
1). 

The overall time of the test 1s limited by the longer T x ; 1f 

is t=5T 1 (where T a is the longer T a value) and the allowed 
observation time 1s T obs , the maximum number of accumulations- 
possible Is N=T obs /(t w +t) T obs /(t w +5T X ). Signal /noise ratio 
is proportional to /N, so that It 1s important to set the pulse 
repetition time x so that t 1s not significantly longer 

20 than 5T 1 ; otherwise information is lost. As mentioned 
previously, a restriction on the separation, Tf, between the two 
pulses 1s that T f should exceed, say, 2T 2 * or 3T 2 *, in order to 
allow the f.i.d. from the first pulse to decay substantially 
before the second pulse is applied. Since T 2 * may vary, for 

25 example, with temperature, it may be important to adjust x f 1n 
the manner described previously to compensate for this. 

In an alternative form of this first embodiment, the pulse 
repetition time t Is made less than 5T 1 (where T x Is now the 
shorter T x value) and the pulse width and/or rf power adjusted to 

30 produce flip angles which are less than and which allow for 
the variation of x/T x with temperature, as described earlier In 
relation to Figure 1. That Is, the excitation is such as would 
produce equal signal strengths at the temperatures at which f x 
and f 2 are the resonant frequencies, the signals are weaker than 

35 when x/T x = 5, but more can be accumulated 1n a given time and a 
lower rf power is required. 
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In a second embodiment of the Invention shown In Figure 3(b), 
case (b), series of excitation pulses are applied at two 
different frequencies f, and : f 2 to excite separate echoes at 
staggered Intervals. The pulse repetition time x can* be equal to 
.5 or less than 5T r (as with the two forms of the first embodiment), 
and. the pulse width and/or rf power may be adjusted if necessary 
to produce flip angles which allow for the variation of t/T, with 
NOR resonance frequency, as described earlier in relation to 
Figure 1. The echoes at x/T t <5 are weaker than when, t/T ^S, but 
10 more can be accumulated In a given time. In order to compensate 
for the variation of T 2 with temperature. It may be important to 
vary the pulse spaclngs between the excitation pulses (x e ) from 
fj to f 2 so that the ratio t e /T 2 would be constant at the 
temperatures at which fj and f 2 are the resonant frequencies. 
15 interleaving of the pulses In the way shown 1n Figure 3(b) is 

possible If t»T 2 . In the "embodiment shown 1n Figure 3(b) the 
pulses are at only two frequencies, f> and f 2 . However, several 
different radlofrequencies could be Interleaved to cover as wide 
a temperature/frequency range as 1s required, consistent with the 
20 band width of the rf probe and receiver. 

In a third embodiment of the Invention, case (c) (not shown), 
a long series of short high-power rf pulses at a single frequency 
and appropriate phases Is used with a spacing x e less than the 
spin-spin relaxation time T 2 ; the responses are then observed as 
25 a train of echoes between each pulse and decay with a time 
constant T 2e , where in favourable cases T 2e >>T 2 . These echoes 
are sampled and accumulated to provide the final signal. The 
whole process is then repeated at one or more different 
frequencies. This embodiment 1s particularly advantageous where 
30 T 2e 1s long. It will be appreciated that any temperature 
variation of the relaxation times, T 1$ T 2 and T 2e , may need to be 
allowed for. 

NQR apparatus for putting the Invention Into effect will now 
be described with reference to Figure 4. Two radio frequency 
35 sources, 1 and 2. provide pulsed rf excitation pulses at 
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frequencies of f x and f 2 , to cover the frequency range produced 
by differing sample temperatures. f x and f 2 are selected 
according to the output of a storage device 3, In which 1s stored 
Information on how the NQR resonance frequency, spin-lattice 
5 relaxation time, T lt spin-spin relaxation times, T 2 , T 2e and 
f.l.d. relaxation time, T 2 * of the sample vary with temperature 
and/or pressure. More frequencies may be necessary 1f the 
frequency range 1s larger than a few tens of kHz. Alternatively, 
a switched frequency synthesizer could be used, provided that the 
10 switching time was much less than the f.l.d. time T 2 * (for RDX, 

" — th1s-quant1ty-1s-about-0r8-ms-at-298K-)-; ; 

The rf sources 1 and 2 have normal and phase-shifted outputs 
(typically 0° and 180°) connected through rf gates 3a,b and 4a,b 
which are controlled by master timing circuitry 5, which also 
15 adjusts the width of each rf pulse according to the excitation 
frequency, the pulse repetition time, t, and the spin-lattice * 
relaxation times stored In the storage device. The rf pulse 
sequences pass to an rf power amplifier 6 of constant output over 
the frequency range of the rf sources, and thence to an rf 
20 probe 7 and rf col 1 (or colls) 8 surrounding the sample. The rf 
probe receives Inputs from the timing circuitry 5 so that the 
tuning elements can be adjusted to take account of- differing 
excitation frequencies f 2 , f 2 , etc., which may be necessary for 
h1gh-Q rf colls. The colls produce a uniform field over the 
25 working region of the probe. The signals generated pass to a 
pre-ampl1f1er 9 and rf amplifier 10 of sufficient band width to 
amplify at constant gain all NQR response signal frequencies 
likely to be encountered. The response signals are detected in 
two separate phase-sensitive detectors lla.b, to generate 
30 in-phase and quadrature outputs. 

The output of gates 3a, b, controlling pulses of radio 
frequency f x (channel 1) 1s connected as a reference signal to 
detector 11a and through a 90° phase-shift network 12 to 
detector lib. Outputs from the detectors lla.b are sampled and 
35 digitised by analogue-to-digi tal converters 13a, b under the 



WO 92/17794 



PCT/GB92/00580 



- 20 - 

control of the timing circuitry 5 and then passed to a digital 
signal processor 14 for display In graphical recorder or video 
display lis. 

The output of gates 4a, b controlling pulses" of radio 
5 frequency f 2 (channel 2) Ms connected as a reference signal to 
the detector 11a and through the 90'-phase-sh1ft network 12 to 
the detector lib. Outputs from the detectors 11a, b are sampled 
and digitised, by the analogue-to-digital converters 13a, b and 
stored In a separate memory .from the signals in channel 1, this 
TO selection being controlled by trigger pulses from the timing 
circuitry 5. After processing 1n the digital signal 
processor 14, both, signals are either separately displayed, or 
* summed before display to provide the output, of the Instrument. 
An alarm (not shown) may then sound If this signal exceeds a 
IS pre-set threshold. 

It will of course be understood that the present Invention 
has been described above purely by way of example, and 
modifications of detail can be made within the scope of the 
Invention. 



WO 92/17794 



PCI7GB92/00580 



- 21 - 
CLAIMS 

1. A method of NQR testing a sample having at least one NQR 
property which varies with a given environmental parameter, 
comprising applying excitation to the sample to excite NQR 

5 resonance and detecting the NQR response signal, the excitation 
being such as to excite an NQR response signal which Is 
detectable during the detection step over a predetermined range 
of the environmental parameter, 

2. A method according to Claim 1 wherein the NQR resonance 
10 frequency of the sample varies with temperature and the 

^ exc1tatton"~ts-~such— as^to-a11ow~for~resonance— frequency— shifts 

caused by ± 20°C temperature variations. 

3. A method according to Claim 1 or 2 wherein the NQR resonance 
frequency of the sample varies with the given environmental 

15 parameter and the excitation Is applied at a single excitation 
frequency, the excitation being sufficient to excite a detectable 
NQR resonance response signal over the range of resonance 
frequency corresponding to the range of the environmental 
parameter. 

20 4. A method according to Claim 1 or 2 wherein the NQR resonance 
frequency of the sample varies with the given parameter and the 
excitation 1s at a sufficient number of excitation frequencies 
that every value of resonance frequency 1n a range corresponding 
to the range of the environmental parameter Is within 0.1% of an 

25 excitation frequency. 

5. A method according to Claim 1 or 2 wherein the NQR resonance 
frequency of the sample varies with the given environmental 
parameter and the excitation is applied at a plurality of 
different excitation frequencies so that a detectable NQR 

30 resonance response signal can be excited over the range of 
resonance frequency corresponding to the range of the 
environmental parameter. 

6. A method according to any of the preceding claims wherein a 
relaxation time of the sample varies with the given environmental 

35 parameter, and the excitation is applied as a series of 



WO 92/17794 



PCT/GB92/00580 



- 22 - 

.excitation pulses, at a plurality of different pulse spaclngs 
and/or so as to produce a plurality of different flip angles, so 
that a detectable NQR resonance response signal can be excited 
over the range of the relaxation time corresponding 'to the range 
5 of the environmental parameter. 

7. A method according to Claim 6 wherein the excitation, is 
applied as a repetitive series of excitation pulses at each of a 
plurality of different excitation frequencies, the pulse 
repetition time being the same for each series, the flip angle 

10~ be1n^1fferent~"for each series. 

8. A method according to Claim 6 or 7 wherein the excitation is 
applied as a series of excitation pulses at each of a plurality 
of different excitation frequencies, each series Including a 
plurality of pulses for generating echoes, the spacing between 

15 each of the plurality of pulses being different for each series. 
.9. A method according to any of the preceding claims wherein the 
excitation is applied as a series of excitation pulses at each of 
a plurality of different excitation frequencies and pulses at the 
different frequencies are applied at staggered spaclngs chosen 

20 according to the variation of T 2 * with the environmental 
parameter. 

10. A method according to any of the preceding claims wherein the 
excitation Is such that the NQR response signal/noise ratio would 
be substantially equal at at leasttwo differing values of the 

25 environmental parameter, and preferably it is such that that 
ratio would be substantially constant over the range of the 
environmental parameter. 

11. A method according to Claim 10 wherein the NQR resonance 
frequency of the sample varies with the given environmental 

30 parameter, and the excitation Is applied at a plurality of 
different frequencies, the number and spacing of the frequencies 
being such that the NQR response signal/noise ratio is 
substantially constant over the range of the environmental 
parameter. . 
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12. A method according to Claim 11 wherein the excitation Is 
applied at n-1 equally spaced frequencies, where n Is the nearest 
Integer satisfying the equation 

ng Av 0 /Av off , 

5 where n£2, Av 0 1s half of the value of the resonance 

frequency range corresponding to the range of the environmental 
parameter, and Av 0 ff 1s approximately equal to the precesslonal 
frequency. 

13. A method according to Claim 10, 11 or 12 wherein a relaxation 
10 time of the sample varies with the given parameter, and the 

"excitation— 1s—app11ed—as~a— series— of— excitation— pulses— at— a 

• plurality of different pulse spaclngs and/or so as to produce a 
plurality of different flip angles, such that the NQR response 
signal /noise ratio is substantially constant over the range of 
15 the environmental parameter. 

14. A method according to Claim 13 wherein the NQR resonance 
frequency of the sample also varies with the environmental 
parameter, the excitation pulses are applied at a plurality of 
different excitation frequencies, and the pulse spaclngs and/or 

20 flip angles associated with the excitation frequencies are such 
as to produce an NQR response signal of substantially equal 
signal/noise ratio at the respective values of the environmental 
parameter at Which the NQR resonance frequency equals the 
respective excitation frequency. 

25 15. A method according to any of the preceding claims wherein the 
excitation is applied at different frequencies to excite separate 
free induction decay responses. 

16. A method according to any of the preceding claims wherein the 
excitation is applied at different frequencies to excite separate 

30 echo responses. 

17. A method according to Claim 15 or 16 wherein the responses 
occur at staggered intervals. 

18. A method according to any of the preceding claims wherein the 
excitation Includes at least one excitation pulse which is so 

35 shaped as to produce a response signal whose strength 1s 
approximately constant over a limited resonance frequency range. 
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J9. Apparatus for NQR testing a sample having at least one NQR 
property which varies with a given environmental parameter, 
comprlslng:- 

storage means In which Is stored Information on' how the or 
each NQR property varies over a predetermined range of the 
environmental parameter; 

means for applying excitation to the sample to excite NQR 
resonance; 

means for detecting the NQR response signal; and 



10 



means, responsive to the Information stored In the storage ( 
..means, for controlling the applying means to excite an NQR 
response signal which is detectable by the detection means over a 
predetermined range of the environmental parameter . 

20. Apparatus according to Claim 19 wherein the storage means': 
15 stores information on how the NQR resonance frequency of the 

sample varies with temperature over a ± 20°C temperature. .range, 
and the controlling means allows for resonance frequency shifts 
oyer such a temperature range. 

21. Apparatus according to Claim 19 or 20 wherein the. storage . 
20 means stores information on how the NQR resonance frequency of 

the sample varies with the given environmental parameter, and the 
controlling means 1s arranged to apply the excitation at a 
plurality of different excitation frequencies so that a 
detectable NQR resonance response signal can be excited over the ' 
25 range of resonance frequency corresponding to the range of the 
environmental parameter. 

22. Apparatus according to any of Claims 19 to 21 wherein the 
storage means stores Information on how a relaxation time of the 
sample varies with the given parameter, and the controlling means 
is arranged to apply the excitation as a series of excitation 
pulses at a plurality of different pulse spaclngs and/or so as to 
produce a plurality of different flip angles, so that a 
detectable NQR response signal can be excited over the range of 
relaxation time corresponding to the range of the environmental 
35 parameter. • 



30 
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23. Apparatus according to any of Claims 19 to 22 wherein the 
control means Is arranged to apply the excitation such that the 

>* NQR response signal /noise ratio would be substantially equal at 
at least two differing values of the environmental parameter, and 
5 preferably such that that ratio would be substantially constant 
over the range of the environmental parameter. 

24. A method of NQR testing for the presence of selected nuclei 
1n a sample, comprising: 

deriving an NQR response signal from the sample 1n such a way 
10 that the signal/noise ratio would be substantially equal at at 

-Heast—two~differ1ng-va-lues—of~a— given— environmental-parameter-;— 

and 

producing an alarm signal from the NQR response signal in 
dependence on whether a predetermined threshold of detection has 
15 been exceeded. 

25. A method according to Claim 24 wherein the nuclei have a NQR 
resonance frequency and/or relaxation time (T lf T 2 , f 2fit T 2 *) 
which vary substantially with variation of the environmental 
parameter, and the at least two differing values are values at 

20 which the resonance frequency and/or relaxation time are 
different. 

26. Apparatus for NQR testing for the presence of selected nuclei 
in a sample, comprising:- 

means for deriving an NQR response signal from the sample in 
25 such a way that the signal/noise ratio would be substantially 
equal at at least two differing values of a given environmental 
parameter; and 

means for producing an alarm signal from the NQR response 
signal in dependence on whether a predetermined threshold of 
30 detection has been exceeded. 

27. A method according to any of Claims 1 to 18, for determining 
the spatial distribution of temperature within a sample, wherein 
a relaxation time (T lf T 2 , T 2e , T 2 *>, preferably spin-lattice 
relaxation time, T lf varies with temperature, and such spatial 

35 distribution Is determined according to the variation of the 
relaxation time within the sample. 
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28. A method for detecting the presence within a larger article 
of a specific substance containing quadrupolar atomic ". nuclei , 
wherein the irradiation can allow for resonance frequency shifts 
caused by ±20°C temperature variations. 
5 29. A method according to Claim 28 wherein the irradiation can 
allow for the resonance frequency shifts by Its frequency being 
modulated or varied. 

30. A method for detecting the presence within a larger article 
of a specific substance containing quadrupolar atomic nuclei, 

TO where1n~the power 'spectrum of the irradiation pulses provides 
substantial power within about 0.1% of any frequency to which NQR 
resonance may be shifted , by any environmental conditions likely 
* to apply to the article. 

31. A method of NQR testing a sample substantially as herein 
15 described with reference to the accompanying drawings. 

32. Apparatus for NQR testing a sample substantially as herein 
described with reference to the accompanying drawings. 

20 PCT 
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